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1^ . ABSTRACT 

The accurate determination of galaxy relative distances is extremely important for the 
empirical calibration of the uncertain metallicity dependence of some standard candles 

■ like Cepheid stars, or for studying the galaxy space distribution and peculiar velocities. 
\J3 ' Here we have investigated the reliability of the widely used /-band Tip of the Red 

■ Giant Branch (TRGB) relative distances for a sample of Local Group galaxies with 

complex Star Formation Histories (SFRs) and Age Metallicity Relationships (AMRs) 

(N ; namely the LMC, SMC and LGS3. The use of the If -band is also discussed. By 

employing theoretical stellar population synthesis techniques, we find that using actual 
determinations of SFR and AMR of the LMC and SMC, their RGB is populated by 
stars much younger (by ~9 Gyr) than the Galactic globular cluster counterparts, on 
which the /-band (and A-band) TRGB absolute magnitude is calibrated. This age 
Qh difference induces a bias in both the photometric metallicity estimates based on the 

comparison of RGB colours with globular cluster ones, and the TRGB distances. The 

■ extent of the distance bias - which is not influenced by the actual value of the TRGB 
I absolute magnitude zero point - is strongly dependent on the specific TRGB technique 

■ applied, and on the assumed 7-band bolometric correction {BCj) scale adopted; the 

correction to apply to the SMC-LMC distance modulus ranges from mag up to 
+0.10 mag. LGS3 is an example of galaxy populated mainly by old stars, so that 
photometric metallicity and distance estimates using globular cluster calibrations are 

' reliable. However, the relative distance moduli between Magellanic Clouds and LGS3 

are affected by the population effects discussed for the LMC and SMC. The correction 
to apply to the LGS3-LMC distance modulus ranges between —0.05 mag and to +0.14 
mag, whereas in case of the LGS3-SMC distance modulus it goes from —0.07 mag to 
+0.04 mag. In case of all three relative distances discussed before, the correction to 
apply to the A-band TRGB distances are larger than the /-band case. 

Our results clearly show that the presence of a well developed RGB in the Colour 
Magnitude Diagram of a stellar system with a complex SFR does not guarantee that 
it is populated by globular cluster-like red giants, and therefore the TRGB method 
for distance determination has to be applied with caution. A definitive assessment of 
the appropriate corrections for population effects on TRGB distances has however to 
wait for a substantial reduction in the uncertainties on the BCi scale for cold stars. 

Key words: Local Group - galaxies: stellar content - galaxies: distances and redshifts 
- Magellanic Clouds - stars: distances 
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1 INTRODUCTION 

The red giant branch (RGB) is one of the most promi- 
nent and well populated features in Colour-Magnitude- 
Diagram (CMD) of stellar populations with ages larger than 
~ 1.5 — 2.0 Gyr. The upper boundary of the RGB locus is 
the so called tip of the RGB (TRGB), populated by stars 
whose electron degenerate core is to start the helium flash. 



At the helium flash the electron degeneracy is lifted and 
quiescent helium burning will proceed in a non degenerate 
core. The helium flash phase is extremely short, and after 
the degeneracy is lifted, stars are located at much lower lu- 
minosities during the quiescent central helium burning. The 
subsequent asymptotic giant branch phase - that basically 
overlaps with the upper RGB and then reaches higher lumi- 
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nosities - is much faster than the RGB one, and therefore 
the TRGB manifests itself as a discontinuity in the differen- 
tial luminosity function of bright red stars in the CMDs of 
intermediate-age to old stellar populations. 

As for the TRGB absolute magnitude, theoretical and 
semiempirical works (see, e.g., Da Costa & Armandroff 1990, 
LFM93, Salaris & Cassisi 1997, 1998) all agree that observa- 
tions of the TRGB in the /-band minimize its dependence 
on age and chemical composition, so that the M™ GB is 
considered to be a very good standard candle for stellar 
populations displaying a RGB in their CMDs. LFM93 pre- 
sented a comprehensive calibration of the TRGB method to 
determine galaxy distances, using results from Da Costa & 
Armandroff (1990 - hereafter DA90). The TRGB has been 
employed to derive distances to objects in the Local Group 
and beyond, out to the Leo I group and Virgo (e.g. Lee, 
Freedman & Madore 1993 - hereinafter LFM93, Sakai et 
al. 1997, Harris et al. 1998, Cioni et al. 2000, Jerjen & Re- 
jkuba 2001, Mendez et al. 2002, Karachentsev et al. 2003, 
McConnachie et al. 2004). Such distance estimates are very 
important for calibrating and testing other distance indica- 
tors, for measuring the local velocity fields and also to deter- 
mine star formation histories of the galaxies under scrutiny. 

Due to the popularity of the TRGB as standard can- 
dle, a series of works have been devoted to improve various 
aspects of the TRGB method. Different techniques for the 
TRGB detection have been put forward, often to improve its 
detection in case of sparsely populated CMDs, or in presence 
of high photometric errors. The most used method is the one 
presented by LFM93, which employs an edge detection algo- 
rithm (zero sum Sobel kernel [-2, 0, +2]) to detect a disconti- 
nuity in the binned luminosity function of bright RGB stars. 
A refinement of this method was presented in Sakai, Madore 
& Freedman (1996), who applied a continuous edge detec- 
tion function to the Gaussian smoothed luminosity func- 
tion. Alternative techniques have been discussed by Cioni 
et al. (2000), Sarajedini et al. (2002), Mendez et al. (2002), 
Frayn & Gilmore (2003), McConnachie et al. (2004). 

On the other hand, various theoretical, semiempirical 
or empirical calibrations of M™ GB have been published by 
Salaris & Cassisi (1997, 1998), DA90, Bellazzini et al. (2001, 
2004). Calibrations of the absolute magnitude of the TRGB 
in infrared bands (J, H, K) are also available (Ferraro et 
al. 2000, Bellazzini et al. 2004), which largely avoid uncer- 
tainties related to the effect of extinction but, as we will see 
later, have a stronger dependence on the metallicity, and 
hence are likely more sensitive to population effects. 

Once a given TRGB calibration is employed - the zero 
points of the existing calibrations span a range of ~0.2 mag 
- relative distance moduli obtained from the TRGB are con- 
sidered to be extremely precise, sometimes with error bars 
of the order of only 0.05 mag or less (e.g. McConnachie et 
al. 2004). This precision has prompted the use of the TRGB 
as a reference standard candle on which to calibrate empiri- 
cally the metallicity dependence of other widely used stellar 
distance indicators, like Red Clump stars (i.e. Pietrzyhski et 
al. 2003) or Cepheid variables (i.e., Sakai et al. 2004). 

A fundamental point to be stressed is that all calibra- 
tions of MJ rgb (or the infrared counterparts) applied to de- 
termine TRGB galaxy distances assume that the observed 
stellar populations are 'old', with ages of the order of the 
Galactic globular cluster ones (12-14 Gyr). This stems from 



the implicit assumption that a well developed RGB is pop- 
ulated by old, globular cluster-like stars. This is possibly a 
very good assumption when observing, e.g., haloes of spiral 
galaxies, but it is not justified in case of composite stel- 
lar populations. Barker, Sarajedini & Harris (2004) have 
recently studied the case of composite stellar populations 
with a fraction of stars with ages of the order of 1.5 Gyr, 
corresponding to the RGB phase transition (Sweigart, Greg- 
gio & Renzini 1990), i.e., when the He ignition happens in 
a progressively less degenerate core. During this transition 
the TRGB moves at lower luminosities with respect to older 
populations and Barker et al. (2004) found that if more than 
30% of the total number of stars created in a given stellar 
population have an age of this order, the distance obtained 
applying the standard TRGB calibrations is overestimated 
by ~ 10 — 20%. However, the authors conclude that the 
galaxy distances obtained from the TRGB method up to 
now have not been affected by this effect, since none of the 
observed stellar populations appears to have a star forma- 
tion history of this kind. 

In this paper we complement the analysis by Barker et 
al. (2004) by studying the effect of the star formation his- 
tory on the TRGB distances, considering specific examples 
of galaxies in the Local Group, e.g. LMC, SMC, and LGS3; 
our results are independent of both the TRGB zero point 
calibration, and the adopted technique to detect the TRGB 
from the stellar luminosity function. We will show how even 
without the presence of a substantial population at the RGB 
phase transition, systematic errors of 0.1-0.2 mag in TRGB 
relative distances might be present when the stellar popula- 
tion under scrutiny is not uniformly as old as the Galactic 
globular clusters. The paper is organized as follows. Sec- 
tion 2 discusses briefly the TRGB method and its calibra- 
tion in both I and infrared passbands, while Sect. 3 presents 
a detailed reanalysis of the TRGB distances to LMC, SMC 
and LGS3 and discusses the biases in their relative distances 
if one follows the traditional use of the TRGB method. Con- 
clusions follow in Sect. 4. 



2 TRGB BRIGHTNESS AS A FUNCTION OF 
AGE AND METALLICITY 

The bolometric luminosity of TRGB stars as derived from 
stellar evolution models is largely determined by the size of 
the He core at the He flash, once the initial chemical compo- 
sition is fixed 1 . Since low mass stars ignite He all with similar 
core masses (slightly increasing for decreasing stellar mass) , 
i\f ™ GB changes by a few hundredths of magnitudes in the 
typical age range of Galactic globular clusters (ages between 
~ 8 and ~ 13 Gyr), as shown in Fig.Qusing the scaled-solar 
models by Girardi et al. (2000). When approaching the RGB 
phase transition M™ GB increases sharply, due to the lift- 
ing of electron degeneracy in the He core, that causes He 
ignition to occur at significantly smaller core masses. 

On the other hand, once the age is fixed, M™ GB de- 
creases for increasing metallicity, in spite of the decrease of 

1 At a fixed initial composition, also presents a very 

small dependence on the envelope composition changes caused 
by the first dredge-up. This dependence however is much smaller 
than the one on the core mass. 
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2.1 TRGB brightness in the /C-band 
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Figure 1. Afbol of the TRGB as a function of age and [Fe/H] 
from the theoretical models by Girardi et al. (2000) 



the He core mass at the TRGB. As shown, e.g., in Kippen- 
hahn & Weigert (1990) on the basis of homology relations, 
this property stems from the fact that at a given core mass 
an increase of the molecular weight in the H-burning shell 
causes an increase of the RGB star luminosity. The net ef- 
fect is a decrease of the TRGB brightness with decreasing 
[Fe/H]; in general, M™ GB oc -0.20 [Fe/H] at a fixed age 
far from the RGB phase transition (say > 3 Gyr) and for 
metallicities well below the solar value. 

DA90 have provided a much used semiempirical calibra- 
tion of M™ GB as a function of [Fe/H] in the Galactic glob- 
ular cluster regime (old ages and [Fe/H] at most ~ —0.7). 
They employed empirical determinations of the bolometric 
magnitude of RGB stars in a sample of clusters of known 
metallicity, assumed a globular cluster distance scale and 
fixed the slope of the Af™ GB vs. [Fe/H] relationship to the 



value given by Sweigart 
obtaining 



Gross (1978) theoretical models, 



Af h J 



-0.19 [Fe/H] -3.81 



(1) 



A previous semiempirical determination by Frogel et 
al. (1983), as well as more recent theoretical and semiem- 
pirical evaluations of this relationship (e.g. Salaris & Cas- 
sisi 1997, 1998, Ferraro et al. 2000) including the results 
displayed in Fig. have a slope very similar to the value 
of the DA90 calibration, although the zero point can differ 
by 0.1-0.2 mag (see, e.g., the discussion in Salaris, Cassisi 
& Weiss 2002). 

Therefore, there is a good general agreement on the 
way M™ GB behaves as a function of age or metallicity. But 
distance determinations using the RGB tip are mostly based 
on I- and Tf-band observations. In the following we will 
separately discuss the calibrations in these two bands. 



Figure [5] displays the run of Mj 



as a function of age 



for various [Fe/H] values. Here we have used as a reference 
the theoretical models by Girardi et al. (2000) , transformed 
to the Tf-band magnitude using the theoretical bolometric 
corrections (BCk) used by Girardi et al. (2000 - dashed 
line), the empirical BCk by Montegriffo et al. (1998 - dot- 
ted line), and the combination of theoretical plus empirical 
BCk (empirical data for effective temperatures T e g <3500 
K) presented in Girardi et al. (2002 - solid line). 

It is easy to notice that in the isT-band the effect of 
metallicity is enhanced with respect to the case of bolometric 
magnitudes. The general behaviour is the same irrespective 
of the use of empirical or theoretical bolometric corrections 
that produces mainly zero point offsets. 

In case of Galactic globular clusters one finds 
semiempirically (e.g. Ferraro et al. 2000) that 
AM K RGB /A[Fe/H] oc~ -0.60, i.e., the metallicity 
dependence is a factor of ~3 higher than for the bolometric 
magnitudes, in agreement with the model predictions. 
This difference stems from the behaviour of BCk, that is 
practically independent of [Fe/H], and strongly affected by 
T e ff (see, e.g., Frogel, Persson & Cohen 1981, Montegriffo 
et al. 1998) in the sense that BCk decreases strongly for 
increasing T e g. This produces fainter magnitudes for the 
hotter TRGBs, which are also the more metal poor and 
fainter ones in A'/boi, hence the increased dependence of 
m trgb Qn [Fe/H] when compared to M™ 03 . 

One can also notice that, due to the behaviour of BCk, 
there is a non negligible dependence of Af™ GB on the age, 
in the sense that a lower age at a given [Fe/H] mimics a 
lower metallicity at a given age. 

In spite of the appeal due to the negligible dependence 



on reddening, M 



TRGB 



is in principle prone to serious un- 



certainties when metallicity and age of the parent stellar 
population are uncertain. 



2.2 TRGB brightness in the J-band 

As shown by e.g., Bellazzini et al. (2004) for Galactic glob- 
ular clusters, the TRGB magnitude dependence on [Fe/H] 
decreases when moving from the 7^-band to the J-band, 
i.e. dM]j RGB /d[Fe/H] oc -0.54 and dAf J RGB /d[Fe/#] oc 
—0.26. In the /-band the trend of the bolometric cor- 
rection BCi with [Fe/H] tends to cancel the dependence 
Fe/H]. As an order of magnitude esti- 
0.19 and recall that 
m trgb = M™ GB - BCi. DA90 found empirically that 
dBCi/d(V - I) oc -0.243 for RGB stars in a sample of 
Galactic globular clusters of various metallicity; Bellazz- 
ini, Ferraro & Pancino (2001) showed empirically that the 
RGB colour in a sample of globulars scales with [Fe/H] 
as d(V - 7)/d[Fe/H] oc 1.162 [Fe/H] + 2.472. This means 
that about [Fe/H]= -2.0 one has dBC//d[Fe/H] oc -0.04 
and about [Fe/H]= -0.7 one gets dBC//d[Fe/H] oc -0.40. 



of M™ GB on . 

mate we use dM^ GB /d[Fe/H] oc 



Recalling the relationship between M™ GB and M™ GB 
we get dMj nGB /d[Fe/H] cc -0.15 about [Fe/H]= -2.0 
and dAf / TRGB /d[Fe/Ji'] oc +0.20 about [Fe/H]= -0.7, with 
an almost negligible dependence on [Fe/H] at intermediate 
metallicities. 



This kind of relationship between M 1 



and [Fe/H] 
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Figure 2. Mk of the TRGB as a function of age and [Fe/H], 
The [Fe/H] values are listed in order of decreasing TRGB A'-band 
brightness. The theoretical models are by Girardi et al. (2000), 
transformed to the i\ -band magnitude using the theoretical bolo- 
metric corrections {BCk) by Girardi et al. (2000 - dashed line), 
the empirical BCk by Montegriffo et al. (1998 - dotted line), 
and the combination of theoretical plus empirical BCk (empir- 
ical data for effective temperatures T a g <3500 K) presented in 
Girardi et al. (2002 - solid line) 
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Figure 3. Mj of the TRGB as a function of [Fe/H], for a set of 
12.5 Gyr isochrones. 



is matched by theoretical models (see e.g., Salaris & Cas- 
sisi 1998) when considering ages typical of Galactic globular 
clusters, on which this calibration is based. To illustrate fur- 
ther this point, Fig.[3]displays the run of Mj RGB as a func- 
tion of [Fe/H] for a typical globular cluster age of 12.5 Gyr, 
using the theoretical models by Girardi et al. (2000) and 
the bolometric corrections by Girardi et al. (2002). A com- 
parison with Fig. [3] immediately highlights the advantage of 
using the /-band instead of K for old stellar populations: 
in fact, for [Fe/H] between -0.7 and -2.0 Mj RGB varies less 
than 0.2 mag, as compared to the ~ 1 mag variation of 

M TRGB 

The preferred TRGB distance methods are based there- 
fore on /-band observations. A very widely used technique is 
the one presented in LFM93, based on Eq. Q supplemented 
by the DA90 empirical relationship 

BCi = -0.243 (V - I) + 0.881 (2) 

and the following relationship between the RGB colour at 
Mr = -3.5 (denoted as (V - /)o,-3.s) and [Fe/H] 

[Fe/H] = -12.64 + 12.6 (V-J) ,-3.b - 3.3 (V-/)S,-3.5- (3) 

The LFM93 method works as follows. The apparent I- 
band magnitude of the TRGB is determined from the po- 
sition of the luminosity function discontinuity, and after an 
estimate of the stellar population reddening, one determines 
the intrinsic colours of the observed RGB stars. A prelimi- 
nary distance modulus is then fixed and employed to deter- 
mine a preliminary metallicity by measuring (V — /)o,-3.5 
and using Eq. [3] With this preliminar [Fe/H] estimate, one 
applies Eq.Qto determine Af™ GB , and Eq.|!2]to determine 
BCi from the intrinsic mean colour of TRGB stars. With 
these two quantities one obtains a second approximation to 
the real distance modulus from 

(m - M) = I ,trgb + Bd - Af b T ol GB . (4) 

This procedure is then iterated until convergency is ob- 
tained. Since the almost vertical nature of the RGB, and 
the weak overall dependency of M™ GB -_BC7 (as discussed 
in case of Galactic globulars on which this calibration re- 
lies) on [Fe/H], convergency is achieved usually within 2-3 
iterations. 

Alternative approaches to TRGB distance determina- 
tions either adopt theoretical BCi and M™ GB values from 
model atmosphere and stellar evolution computations to de- 
termine directly a theoretical relationship A/? RGB -[Fe/H] 
(as in Salaris & Cassisi 1998) for globular cluster-like ages, or 
derive empirical M/ RGB -[Fe/H] relationships by detecting 
the TRGB in Galactic globular clusters with (empirically) 
known distances and metallicities (e.g. Bellazzini, Ferraro & 
Pancino 2001, Bellazzini et al. 2004). 



3 APPLICATION TO COMPOSITE STELLAR 
POPULATIONS 

As discussed extensively in the previous section, all TRGB 
absolute magnitude calibrations applied for determining ex- 
tragalactic distances are based on Galactic globular cluster 
stars (observations, models or both). Any difference between 
the calibrating stellar populations and the observed ones 
may have an impact on the inferred TRGB distances. As we 
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will demonstrate shortly, the presence of a well populated 
RGB in a generic CMD does not automatically imply that 
the stellar population under scrutiny is equivalent to Galac- 
tic globular cluster populations. In case of spiral galaxies, ob- 
serving far in the halo will very likely minimize population 
differences, but when working with, e.g., irregular galaxies 
(like many of the components of the Local Group) there is 
not a clear-cut distinction between the homogeneously old 
stellar component and younger ones. 

In the following we will discuss these population effects 
on the TRGB distances to the LMC and SMC by using 
population synthesis simulations of the galaxy stellar popu- 
lations, in the same vein as in Girardi & Salaris (2001). 

3.1 LMC 

3.1.1 The simulated CMD 

We have computed synthetic CMDs of the LMC stellar pop- 
ulations using the Star Formation History (SFH) and Age 
Metallicity Relation (AMR) displayed in Fig. 0] which are 
typical of the LMC bar stellar populations (see discussion 
in Girardi & Salaris 2001, Salaris, Percival & Girardi 2003). 
The synthetic CMDs have been computed using the same 
code employed by Girardi & Salaris (2001), and we made 
use of the stellar models by Girardi et al. (2000) and colour 
transformations by Girardi et al. (2002). We have included in 
the simulation a Gaussian la photometric error of 0.02 mag 
(typical error at the magnitudes of LMC and SMC TRGB, 
see, e.g., the discussion in Barker et al. 2004), and a Gaussian 
la spread of 0.1 dex around the [Fe/H] values displayed in 
fig. 0] but the following results are completely independent 
of the photometric errors and [Fe/H] dispersion. 

The top panel of Fig. shows the Mi — (V — I) CMD 
of the synthetic LMC RGB population. It is evident, in this 
figure, the position of the TRGB at about Mi ~ —4. 

3.1.2 The apparent discrepancy in mean [Fe/H] 

Overplotted in Fig. are the RGBs of two globular cluster 
isochrones (age of 12.5 Gyr) with the labelled [Fe/H] of —0.9 
and —1.5. They delimit the bulk of RGB stars found in the 
simulation. By comparing the (V — I) colours of the globular 
cluster isochrones and the LMC population, we can derive, 
in a way similar to the use of Eq. |3] the mean [Fe/H] value 
for the RGB stars in the LMC. 

A first fundamental point has to be noticed at this stage. 
Our synthetic LMC population has been computed using the 
scaled solar models by Girardi et al. (2000); the LMC chemi- 
cal evolution models by Pagel & Tautvaisiene (1998) and the 
spectroscopic observations by Hill et al. (2000) and Smith et 
al. (2002) show that at low metallicity the a-elements (i.e. 
O, Ne, Mg, Si, S, Ca, Ti) are mildly enhanced with respect to 
Fe, i.e., [a/Fe]~0. 1-0.2. However, the relevant [Fe/H] range 
of the RGB stars will turn out to be - as thoroughly dis- 
cussed in the following - such that [a/Fe]~0 according to the 
results mentioned before, and our use of scaled solar models 
for the simulation of LMC red giants is fully justified. 

The Galactic globular cluster RGB isochrones employed 
to determine the photometric metallicity of our LMC syn- 
thetic population, following the LFM93 method, must how- 
ever take into account the fact that [a/Fe]~0.3 in globu- 
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Figure 4. SFR and AMR adopted for our LMC simulation. The 
SFR (bottom panel) is the one derived by Holtzman et al. (1999, 
his figure 4) from HST photometry of a field in the LMC bar. 
The AMR (top panel) comes from the bursting chemical evolution 
model of Pagel & Tautvaisiene (1998), and agrees well with AMRs 
derived from observations of LMC star clusters (see e.g. Hill et 
al. 2000). 



lar cluster stars (see, e.g., Salaris, Chieffi & Straniero 1993, 
Carney 1996 and references therein). For metallicities typi- 
cal of the Galactic halo a-enhanced theoretical stellar mod- 
els are well reproduced by scaled solar ones with the same 
total metallicity [M/H] (Salaris et al. 1993), and there- 
fore the globular cluster isochrones displayed in Fig. 
are the scaled solar ones by Girardi et al. (2000) with 
the appropriate choice of the total metallicity that mim- 
ics [a/Fe]=0.3, i.e. [M/H] = [Fe/H]+0.2 (for a scaled solar 
mixture [M/H] = [Fe/H]). 

The mean [Fe/H] obtained comparing the mean colour 
of our synthetic LMC population at Mi — —3.5 (the ref- 
erence magnitude used in the LFM93 method) with the 
colours of individual globular cluster isochrones, of various 
[Fe/H] and at the same Mi, is [Fe/H]=-1.08. This is in very 
good agreement with values obtained by LFM93 and Sakai 
et al. (2000) from VI LMC observations. We also redeter- 
mined the mean [Fe/H] of LMC giants employing photome- 
try of the LMC OGLE field 6 (Udalski et al. 2000) and the 
method by LFM93. After correcting for the effect of extinc- 
tion using the reddening maps by Udalski et al. (1999) we 
obtained a mean [Fe/H] = — 1.15. 

This [Fe/H] value obtained from the synthetic CMD 
and in agreement with estimates from observed CMDs is 
however greatly discrepant from the results of spectroscopic 
observations of bright RGB stars in the LMC bar published 
by Cole et al. (2000), who found a mean [Fe/H] = -0.60. It is 
important to notice that if we apply the LFM93 technique 
and calibrations, this [Fe/H] difference, taken at face value, 
would produce a systematic change in the distance modulus 
by ~0.10 mag, because BCi in Eq.0]is fixed by the observed 
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Figure 5. Simulation of the LMC bar population using the SFR 
and AMR of Fig. HI The top panel shows the simulated Mj vs. V — 
I CMD detailing the upper part of the RGB (dots), and compares 
it with two globular cluster isochrones for [Fe/H] equal to —1.5 
and —0.9, respectively (continuous thick lines). The bottom panel 
shows the metallicity distribution of the simulated upper RGB 
stars. 



Figure 6. Luminosity functions for all red stars [V — I > 1) in 
the LMC simulation (heavy line), and separated in four age bins 
(light lines). Notice that the TRGB discontinuity, detectable at 
about Af(, ; = —3.6, is mainly caused by stars in the 2 — 4 and 
4 — 6 age bins. 



(V - J) colour, and M™ GB changes by 0.10 mag for a 0.5 
dex change of [Fe/H]. 

The bottom panel of Fig. |S] shows the [Fe/H] distri- 
bution for the stars in our simulation that populates the 
upper RGB. Here we find a mean [Fe/H] =— 0.61, in excel- 
lent agreement with the spectroscopic result of —0.60. This 
means that also the synthetic LMC population displays the 
~0.5 dex discrepancy between the RGB [Fe/H] inferred from 
the RGB colour and the true value (the spectroscopic one 
in case of real LMC stars). 

The reason for this discrepancy is twofold. The first 
cause is related to the age distribution of LMC stars com- 
pared to the age of globular clusters. Figure HJ displays the 
binned luminosity function (bins 0.04 mag wide) of all bright 
red stars (V — I >1.0) in the LMC 2 and the partial lumi- 
nosity functions for different age ranges. The discontinuity 
corresponding to the TRGB is clearly visible, and for what 
follows the technique used to identify the location of the 
discontinuity is completely irrelevant. 

The mean age of bright RGB stars is ~4 Gyr, much 
lower than the age (12-13 Gyr, see e.g. Salaris & Weiss 2002) 
of the Galactic globular clusters on which the RGB colour- 
[Fe/H] relationship is calibrated. Although the RGB colour 
is weakly sensitive to the age of the parent stellar population, 
an age difference of ^9 Gyr has an appreciable influence, in 
the sense that the younger LMC population is bluer than 

2 We display the luminosity function in the bolometric magni- 
tude since the LFM93 method employs a theoretical calibration 
of MTJ$ GB as the standard candle of choice. Identical results are 
obtained from the luminosity function in either the I- or K-haad. 



the globular cluster one at fixed chemical composition. If 
we shift the mean colour of the LMC RGB to the red, to 
account for the effect of a 9 Gyr age difference with respect 
to a typical globular clusters RGB, we would determine a 
photometric [Fe/H] value that is ~0.3 dex lower than the 
intrinsic one. 

After the correction for the age effect there are still 0.2 
dex of mismatch in [Fe/H] to be accounted for. This remain- 
ing difference is related to the different metal mixture of the 
calibrating globulars. Since a given [Fe/H] value coupled to 
[a/Fe]>0 corresponds to a higher [M/H] than for a scaled 
solar mixture, RGB stars in the LMC appear to be bluer 
than the counterpart in globulars with the same [Fe/H]. 
This causes an underestimate of the LMC stars [Fe/H] of 
about 0.2 dex if [a/Fe]~0.3 for Galactic globular clusters, 
and in this way the full difference between photometric and 
real [Fe/H] is explained. 

Let us for a moment suppose that the underestimate of 
[Fe/H] for the LMC RGB is due only to this latter effect, and 
denote with [Fe/H]' the value obtained from Eq. ^increased 
by 0.2 dex. We know that M™ GB is determined by the total 
stellar metallicity [M/H] (e.g. Salaris et al. 1993), and that 
the calibration given by Eq.^is based on a-enhanced glob- 
ular cluster stars, for which [M/H] = [Fe/H]4-0.2. This means 
that the use of the photometric [Fe/H] given by Eq. in- 
stead of [Fe/H]' does provide the correct M™ GB for the 
LMC stars, because at a given [Fe/H] their [M/H] is lower 
by 0.2 dex than the globular cluster counterpart, and this 
mismatch is compensated for by the use of the 'wrong' pho- 
tometric [Fe/H]. The bottom line is that, from the point 
of view of the expected absolute magnitude of the TRGB 
for LMC stars, the 0.2 dex [Fe/H] mismatch due to the a- 
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Figure 7. Difference between the theoret- 
ical M™ GB for t=12.5 Gyr, scaled so- 
lar [Fe/H]=— 1.38 (corresponding to [Fe/H]~ 
— 1.6 for the typical Galactic globular clus- 
ter (^-enhancement) and the theoretical val- 
ues predicted for the ages and scaled solar 
[Fe/H] values displayed. The underlying the- 
oretical models are from Girardi et al. (2000). 
Panel a), b), c) and d) show, respectively, the 
results using the Yale transformations, West- 
era et al. (2002) transformations, Girardi ct 
al. (2002) transformations, and the transfor- 
mations used in Girardi et al. (2000). 



enhancement in Galactic globular clusters does not play a 
role in the LFM93 method. 



3.1.3 Implications for the TRGB distance to the LMC 
using the LFM93 method 

A few facts can be derived from the above discussion: 

(1) The TRGB discontinuity (in Mboi, but the same is 
true for the I- and if-band) in the LMC is determined by 
stars with ages of the order of ~4 Gyr. Older and younger 
RGB stars, although certainly present, are either too few 
(in the case of old stars), or do not contribute significantly 
to the detected TRGB discontinuity (in the case of stars 
younger than 2 Gyr), or both. 

(2) Since most of the RGB stars in the LMC are 
intermediate-age, it is wrong to determine their [Fe/H] by 
simply comparing their V — I colours with those of globular 
cluster giants. 

(3) When using the LFM93 method in the LMC, a num- 
ber of additional steps have to be followed: Whereas BCi is 
fixed by the observations, M™ GB has to be corrected for a 
0.3 dex metallicity increase effect (that increases the M™ GB 
at fixed age) and for the age decrease with respect to the 
globular cluster ages (this effect is however very small), as 
predicted by the theoretical models. Applying these correc- 
tions, the LMC distance modulus increases by 0.05 mag with 
respect to results with the original LFM93 method and cal- 
ibrations. 

A particular remark should be made regarding item (2). 
As seen in Fig. 4, in the LMC the SFR has been variable 



and presents a gap from 7 to 9 Gyr, that nearly coincides 
with the usual separation betwen populations of intermedi- 
ate (~ 1 to 7 Gyr) and old ages (>9 Gyr). At intermediate 
ages, the mean SFR is just 3.1 times higher than the mean 
one at old ages. It may seem surprising that this moderate 
excess of the SFR at intermediate ages has caused the RGB 
to be dominated by ~4 Gyr old stars, with just a minor 
contribution from old ages (see Fig. 6). Actually, there is no 
surprise in this result: The number of evolved stars of a given 
kind and age is proportional to both the SFR at that age, 
and to their intrinsic lifetime and birth rate. The birth rate 
is a function that strongly decreases with the stellar age, 
and this would cause the predominance of intermediate-age 
stars even for the case of constant SFR. This issue has been 
thoroughly discussed, in the case of red clump stars, by Gi- 
rardi & Salaris (2001); the situation for RGB stars close to 
the TRGB is a very similar one. 

Therefore, the dominance of the intermediate-age RGB 
is not a result specific to the LMC: it is the case for any 
galactic component in which the star formation, at interme- 
diate ages, has been significant compared to old ages. This 
includes the discs of spirals, most of the irregulars, and a 
significant fraction of the dwarf spheroidal galaxies. 



3.1.4 Implications for the TRGB distance to the LMC 
using other methods 

Other methods usually consider a constant M™ GB for 
'metal poor' stellar populations taken from calibrations 
based on Galactic globulars. For example, Frayn & 
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Gilmore (2003), Karechentsev et al. (2003), McConnachie 
et al. (2004), use a constant Af/ RGB value based on globu- 
lar cluster calibrations, and apply it to the observed RGB 
population of the parent galaxy. Of course, in these cases we 
have again the assumption that the RGB population, or at 
least the RGB stars that cause the detected TRGB discon- 
tinuity, are old and metal-poor. As we have just shown, this 
is not the case for the LMC bar population. 

To check the bias in the LMC distance obtained with 
this method one can use theoretical isochrones differentially, 
once we have identified - by means of simulations like the 
one of Fig. |S] - the typical ages and metallicities of TRGB 
stars. 

First, we can directly compare M™ GB obtained from 
the isochrones for a typical halo globular cluster metallic- 
ity [Fe/H]=-1.6 and an age of 12.5 Gyr, with the Mj RGB 
obtained using the true [Fe/H] and true age of the TRGB 
stars in the LMC. The difference between the two Mj RGB 
values, hereafter defined as 



amJ 



Mi 



TRGB (12.5Gyr, -1.6) - M 7 TRGB (trwe) 



(•») 



provides the correction we have to apply to our distance 
determination. We have performed this exercise, based on 
the synthetic data of Fig. |S] 

Somewhat surprisingly, we find that the results depend 
in a non-negligible way on the adopted set of transformations 
between the theoretical (Mb i,T ff) and the observational 
(Mi, V — I) plane. The results can be summarized as follows: 

• Using Girardi et al. (2000) transformations the LMC 
distance has to be increased by 0.07 mag; 

• using Girardi et al. (2002) transformations the LMC 
distance has to be decreased by 0.14 mag; 

• using Yale transformations (Green 1988) the LMC dis- 
tance is unchanged; 

• using Westera et al. (2002) transformations the LMC 
distance is unchanged. 

Figure |7| illustrates the behaviour of the quantity 
AM ™ GB as a function of age and metallicity, for the above- 
mentioned sets of transformations. We recall that they are 
derived from libraries of stellar spectra comprising different 
metallicities. Since studying the reasons for the differences 
among these curves is beyond the scope of the present pa- 
per, we limit ourselves to a few comments. In case of old 
metal poor populations, differences in AM™ GB are limited 
to ~ 0.1 mag, but anyway cannot be considered negligi- 
ble. At high metallicities and old ages, as soon as the T B g 
value of TRGB stars becomes lower than ~ 4000 K, the dis- 
crepancies among the curves become dramatic. In fact, the 
transformations in the range Teff <4000 K are (1) subject to 
larger variations due to the appearance of molecular bands 
in stellar spectra, and (2) dealt with differently by different 
authors. A value T c s — 3500 K is the minimum temperature 
for the ATLAS9 theoretical model atmospheres (see, e.g., 
Kurucz 1993), which constitute the backbone of the trans- 
formations by Westera et al. (2002) and Girardi et al. (2000, 
2002). Below this T e g (and to a certain extent also in the 
range ~ 3500 — 4000 K), the different transformations rely 
on sets of empirical spectra and colour-Toff relations, which 
are more prone to be implemented in different ways. There- 
fore, it is no surprise that AM™ GB heavily depends on the 
adopted transformations in case of high metallicities. 



It is also important to notice that the theoretical uncer- 
tainties in AAf/ RGB cannot be avoided by using the empir- 
ical relationship to compute BCi employed by the LFM93 
technique (see Eq. since it has a non negligible la scat- 
ter of 0.057 mag that DA90 attribute to observational errors. 
Moreover this empirical BCi calibration by DA90 does not 
reach the true mean [Fe/H] of the LMC RGB stars, its up- 
per limit of validity being [Fe/H] = — 0.7 (the metallicity of 
47 Tuc). 

Could one avoid the uncertainties in the transforma- 
tions, by deriving a TRGB distance using infrared instead 
of /-band observations? As indicated by Fig. |2] different 
sets of transformations provide almost the same behaviour 
of M™ GB as a function of age and metallicity, thus sug- 
gesting that AA/™ GB could be derived more reliably than 
am trgb The problem in thig 

case is related to the high 
dependence of A'/£ RGB on age and metallicity: The compu- 
tation of the correction AM™ GB requires the knowledge of 
age and metallicity distribution of RGB stars, any error in 
these quantities are likely to produce substantial errors in 
AAf™ GB . Since the calibration of M™ GB is based on glob- 
ular clusters - located at one extreme of the age-metallicity 
plane allowed to RGB populations - biases are expected to 
be high. The maximum bias is obtained when using the K- 
band magnitude of the TRGB coupled to a semiempirical or 
theoretical Af™ GB calibration based on globular clusters, 
and a calibration (e.g. Ferraro et al. 2000) (V - A")-[Fe/H] 
or (J — A")-[Fe/H] obtained from globulars. By using differ- 
entially the theoretical isochrones in infrared passbands, the 
[Fe/H] and age biases discussed before would produce a bias 
of ~0.2 mag (true distance larger than the value obtained 
from the globular cluster calibration) in the TRGB distance, 
when considering a globular cluster age of 12.5 Gyr. In ad- 
dition, for the if-band - as mentioned before - the precise 
age of the Galactic globulars plays also a role. 



3.2 SMC 

After discussing in detail the case of the LMC, we now ex- 
amine the SMC. Many of the considerations made for the 
LMC apply also to the SMC, and will not be repeated for 
the sake of conciseness. 

For our simulations, we select the 'global' SFR and 
AMR derived by Harris & Zaritski (2004) from UBVI pho- 
tometry of the entire main body of the SMC, and shown in 
Fig- [HI Their SFR is certainly the best available estimate for 
the SMC, whereas the AMR relation is in good agreement 
with those derived from observations of SMC star clusters. 
Scaled solar models are again appropriate for the relevant 
metallicities of the RGB stars, according to the discussion 
in Pagel & Tautvaisiene (1998 - their Fig. 9 and discussion 
in Sect. 3). 

The simulation is presented in Fig. [5] The photometric 
mean [Fe/H] , as derived from comparing the RGB with those 
of old isochrones and considering their a-enhancement, is 
[Fe/H]=-1.4. From OGLE-II photometry of SMC OGLE 
field 3, dereddened, we get a similar value, [Fe/H] = — 1.5. 

The true mean RGB metallicity, as derived from the 
simulation itself, is [Fe/H] =— 0.82. Again, we have a discrep- 
ancy of 0.6 dex between the two values of [Fe/H] , that can be 
explained as follows: 0.2 dex are due to the a-enhancement, 
0.4 dex are due to the mean younger age of SMC RGB stars 
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Figure 8. SFR and AMR adopted for our SMC simulation. Both 
SFR and AMR are taken from Harris & Zaritsky (2004; their 
figures 7 and 12, respectively). The SFR represents the global 
star formation in the SMC. The AMR agrees well with those 
derived from observations of SMC star clusters. 
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Figure 9. The same as Fig.0 but for the SMC, using the AFR 
and AMR of Fig. [S] The overplotted globular cluster isochrone 
(thick continuous line) has [Fe/H]=— 1.4. 



with respect to old globular clusters. In fact, the mean age 
of bright RGB stars in the SMC is of just ~3.7 Gyr. 

The TRGB discontinuity is mainly determined by stars 
with ages ~4 Gyr, and not by the oldest component with 
globular cluster-like ages, as is evident by looking at the 
partial LFs depicted in Fig. 1101 



Figure 10. The same as Fig.0 but for the SMC simulation. The 
age ranges for the older ages are slightly different, to highlight the 
different age distribution of SMC stars with respect to the LMC. 



Using the LFM93 method, the effect of the [Fe/H] dis- 
crepancy is an increase of 0.07 mag for the SMC distance. 

Using the 'constant AfJ RGB ' methods, the results again 
depend on the selected set of transformations. Using Girardi 
et al (2002) ones, the SMC distance has to be decreased by 
0.04 mag. With the other transformations tested in Fig. |7| 
the SMC distance is basically unchanged. In the Tf-band, 
the AM™ GB amounts to ~0.2 mag, for all the different 
sets of transformations. 

Comparing the LMC and SMC cases, we can conclude 
the following: In both cases, the use of LFM93 will lead to a 
small underestimate of their distance moduli, amounting to 
0.05 for the LMC, and 0.07 mag for the SMC. These small 
errors are still lower than the overall accuracy of the TRGB 
methods (defined by the uncertainty in the zero points, of 
0.1 — 0.2 mag), but are not completely negligible. Moreover, 
the errors seem to be systematic in the sense of decreasing 
distances to the galaxies that have had significant star for- 
mation at intermediate-ages. Fortunately enough, there is 
no impact on the relative SMC-LMC distance (but one has 
always to keep in mind that the BCi relationship by DA90 
is used outside its range of validity for the LMC TRGB). 

When using the 'constant M™ GB ' methods, instead, 
the results depend on the set of transformations one uses to 
derive the quantity AMj RGB . Using the Girardi et al. (2002) 
transformations, for instance, one derives that the straight 
use of a constant A4j RGB value lead to errors in distance 
moduli amounting to —0.14 mag for the LMC, and —0.04 
mag for the SMC. This would imply an error of 0.10 mag 
(true distance being longer) in the derived relative SMC- 
LMC distance. This is a significant error, that may bias the 
calibration of other standard candles as a function of the 
metallicity. On the other hand, if one uses Westera et al. 
(2002) transformations to derive AAf7 RGB , both SMC and 
LMC distance moduli remain unchanged, as well as their rel- 
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Figure 11. SFR and AMR adopted for our LGS3 simulation, 
from Miller et al. (2001). 



Figure 12. Age distribution of the RGB stars obtained from the 
LGS3 simulation. 



ative distance and the calibration of other standard candles 
on the Magellanic Clouds. 



amounting to a significant 0.14 mag on the relative distance 
modulus. 



3.3 LGS3 

As an example of galaxy where the globular cluster TRGB 
calibration is appropriate, we have considered the case of 
LGS3. Miller et al. (2001) provide a global SFR and AMR 
for this galaxy, displayed in Fig. 1111 A main burst of star 
formation happened at the beginning of the galaxy life, after 
which the SFR level has been very low and almost constant. 
The spread around the mean AMR is also accounted for, 
following Miller et al. (2001) estimates. 

Fieure tT^l displavs the age distribution of RGB stars ob- 
tained from our simulation. The main component is made 
of the oldest objects formed during the strong initial burst 
of star formation. RGB stars born after the initial burst 
amount only to about half of the oldest objects. Since the age 
and metallicity ([Fe/H]=~ —1.3 for the majority of LGS3 
RGB objects) of the bulk of the LGS3 RGB stars is com- 
parable to the globular cluster counterparts, the standard 
calibrations of the TRGB absolute magnitude can be ap- 
plied with confidence to this galaxy. 

More in detail, when employing the LFM93 method a 
0.2 dex underestimate of [Fe/H] is possible, if LGS3 stars 
formed all with a scaled solar metal distribution (as assumed 
in our simulation). However, as seen before, this does not af- 
fect the TRGB distances with LFM93 technique. Also when 
employing the 'constant M™ GB ' methods, the true distance 
is recovered. 

The relative distance modulus LGS3-LMC is affected at 
the level of 0.05 mag by the population effects discussed be- 
fore when using the LFM93 method. The 'constant M™ GB ' 
methods provide different results, depending on the set of 
bolometric transformations adopted, the maximum effect 



4 CONCLUSIONS 

The previous analysis of the TRGB distances to the galax- 
ies LMC, SMC and LGS3 allowed us to point out a series 
of uncertainties affecting the current TRGB distance deter- 
minations and RGB photometric metallicity estimates for 
composite stellar populations. 

All TRGB distance methods rest on the hypothesis 
that the RGB populations observed in external galaxies are 
as old as the stars in Galactic globular clusters, on which 
the TRGB calibrations are based. This is however not al- 
ways true, as we have demonstrated for the specific cases 
of the LMC and SMC, whose RGB populations have ages 
of the order of 3-4 Gyr, instead of 12-13 Gyr typical of the 
Galactic globulars. This age difference (and eventual dif- 
ferences in the [a/Fe] ratio between the Galactic globular 
clusters and external galaxies) affects the TRGB distances 
and [Fe/H] estimate of the parent RGB stars, as has been 
already pointed out by Davidge (2003), Rizzi et al. (2003), 
Salaris et al. (2003). In case of LMC and SMC we found 
that [Fe/H] estimates based on the comparison of the RGB 
colour with the globular cluster counterparts are too low by 
~0.5 dex due to the effect of both age (~0.3 dex) and dif- 
ferent heavy element distribution (~0.2 dex). These effects 
explain the discrepancy between the photometric metallicity 
of the LMC RGB determined from a comparison with the 
colours of globular cluster RGB stars, and the spectroscopic 
measurements . 

The younger age of a galaxy RGB population with re- 
spect to globular cluster ages may induce a bias also in the 
/-band TRGB distance, whose extent is strongly dependent 
on the adopted BCi scale. TRGB relative distances between 
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Magellanic Clouds-type galaxies and galaxies like LGS3 pop- 
ulated by mainly globular cluster-like stars can be apprecia- 
bly affected by these population effects, up to ~0.15 mag, 
depending on the BCi scale adopted. More in detail, the cor- 
rection to apply to the LGS3-LMC distance modulus ranges 
between —0.05 mag and +0.14 mag, whereas in case of the 
LGS3-SMC distance modulus it ranges between —0.07 mag 
and +0.04 mag. As for the distance modulus SMC-LMC, 
the correction due to population effect ranges from mag 
up to +0.10 mag. 

The use of K-bwaA TRGB calibrations - the if-band 
bolometric correction scale appears to be more secure as far 
as the trend with T e fr is concerned - does not help much, 
due the strong dependence of M™ GB on both [Fe/H] and 
age; biases of the order of 0.2 dex can be present in K- 
band TRGB distances when using calibrations of the TRGB 
method based on Galactic globular cluster stars. 

In order to isolate the most metal poor and oldest RGB 
component to which apply safely the TRGB globular clus- 
ter calibration, one might try to apply a colour cut to the 
RGB data, i.e. selecting stars bluer than a given (V — I) 
value, as suggested by our referee. We have however verified 
in the specific cases of the SMC and LMC that this proce- 
dure does not achieve in general the desired results. More in 
detail, we have divided the RGB with -3.70 < Mi < -4.05 
into three colour intervals, i.e. 1.3 < (V — I) < 1.5 (blue), 
1.5 < (V — I) < 1.7 (middle) and 1.7 < (V - I) < 1.9 (red) 
and studied the associated distributions of age and [Fe/H] 
values. For SMC stars we find that the oldest population is 
spread between the red and middle interval, whereas for the 
LMC the oldest population is located in the blue interval. 
This difference reflects the AMR of the two galaxies. In case 
of the SMC the AMR used in these simulations is approx- 
imately flat for ages down to ~ 3 Gyr; a sizable fraction 
of RGB stars has therefore a very similar metallicity, and 
the oldest objects are necessarily located towards the red 
side of the RGB. For the LMC the oldest stars are also the 
most metal poor and therefore they tend to cluster at the 
blue side of the RGB. It is evident from this example that 
the same colour cut applied to different galaxies may imply 
very different distributions of ages and [Fe/H] due to dif- 
ferent star formation and chemical enrichment histories; we 
have therefore to conclude that there is no general rule for 
selecting a particular kind of RGB stellar population based 
on colour cuts along the RGB CMD. 

As for the possibility to detect the TRGB of the oldest 
populations in the LMC and SMC using the colour analysis 
described before, the following has to be noticed. The oldest 
SMC population is generally superimposed to the younger 
and more metal rich one, that determines the level of the 
TRGB discontinuity. As for the LMC, the bluest RGB stars 
in a narrow colour range (1.3 < (V — I) < 1.5) can in prin- 
ciple provide a reasonably clean sample of an old globu- 
lar cluster-like population, but their CMD location overlaps 
with a sizable number of bluer objects belonging to differ- 
ent evolutionary phases of much younger ( < 0.5 Gyr) pop- 
ulations, that can be clearly seen in the upper left part of 
Fig. 5. The large fraction of these non-RGB objects popu- 
lating the same colour bin may hamper a clear detection of 
the LMC old metal poor TRGB. 

To conclude, our results show that the presence of a 
well developed RGB in the Colour Magnitude Diagram of 



a stellar system with a complex SFR does not guarantee 
that the RGB is populated by globular cluster-like red gi- 
ants, nor that its TRGB is determined by them; this means 
that the TRGB method for distance determinations has to 
be applied with caution to all galaxies which present signa- 
tures of intermediate- age stars. A definitive assessment of 
the appropriate corrections for population effects on TRGB 
distances has however to wait for a substantial reduction of 
the uncertainties on the BCi scale for cold stars. 
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